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ABSTRACT 

The final stage of terrestrial planet formation consists of the cleanup of residual planetesimals after 
the giant impact phase. Dynamically, a residual planetesimal population is needed to damp the high 
eccentricities and inclinations of the terrestrial planets to circular and coplanar orbits after the giant 
impact stage. Geochemically, highly siderophile element (HSE) abundance patterns inferred for the 
terrestrial planets and the Moon suggest that a total of about 0.01 M e of chondritic material was 
delivered as 'late veneer' by planetesimals to the terrestrial planets after the end of giant impacts. Here 
we combine these two independent lines of evidence for a leftover population of planetesimals and show 
that: (1) a residual population of small planetesimals containing 0.01 M§ is able to damp the high 
eccentricities and inclinations of the terrestrial planets after giant impacts to their observed values. 
(2) At the same time, this planetesimal population can account for the observed relative amounts of 
late veneer added to the Earth, Moon and Mars provided that the majority of the accreted late veneer 
was delivered by small planetesimals with radii < 10 m. These small planetesimal sizes are required to 
ensure efficient damping of the planetesimal's velocity dispersion by mutual collisions, which in turn 
ensures sufficiently low relative velocities between the terrestrial planets and the planetesimals such 
that the planets' accretion cross sections are significantly enhanced by gravitational focusing above 
their geometric values. Specifically we find, in the limit that the relative velocity between the terrestrial 
planets and the planetesimals is significantly less than the terrestrial planets' escape velocities, that 
gravitational focusing yields a mass accretion ratio Earth/Mars ~ (p® / p m ars)(R® / Rmars) 4 ~ 17, 
which agrees well with the mass accretion ratio inferred from HSEs of 12-23. For the Earth-Moon 
system, we find a mass accretion ratio of ~ 200, which, as we show, is consistent with estimates of 
150-700 derived from HSE abundances that include the lunar crust as well as mantle component. We 
conclude that small residual planetesimals containing about ~ 1% of the mass of the Earth could 
provide the dynamical friction needed to relax the terrestrial planet's eccentricities and inclinations 
after giant impacts, and also may have been the dominant source for the late veneer added to Earth, 
Moon and Mars. 

Subject headings: planetary systems: general — planets and satellites: formation — solar system: 
formation 



1. INTRODUCTION 

Terrestrial planet formation is generally consid- 
ered to consist of three main stages. The first 
stage consis ts of the formation of small planetesi- 
mals (e.g. Chia nlTfc Youdinl l2010f h the second stage 
of the coagulation of these small planetesi mals into 
roughly Mars-sized protop lanets (e.g. Ilda fc Makinol 
119931 : IWeidenschilling et al.l[l997l ). and the third stage 
is comprised of collisions of a few dozen protoplanets, 
called giant impacts (e.g. lAgnor et al.l 119991 : iChambersI 
l2001f ). However, several lines of evidence suggest that a 
significant amount of mass was left in planetesimals at 
the end of the giant impact phase in the terrestrial planet 
region. This therefore argues for an additional and final 
stage of terrestrial planet formation, which consists of the 
clean-up of the leftover planetesimals. This final stage of 
terrestrial planet formation is the focus of the present 
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paper. 

Evidence for a significant population of planetesimals 
in the terrestrial planet region at the end of the giant 
impact phase comes from two different areas of research: 

Geochemical evidence suggests that the Earth accreted 
chondritic materials equivalent to about 0. 3%-0.7% of it s 
total mass after the end of giant impacts (Walker 2009). 
The evidence for this 'late veneer' comes from highly 
siderophile elements (HSEs) that are geochemically char- 
acterized as having a strong tendency to partition into 
metal relative to silicates. Hence the silicate portions 
of terrestrial planets with rocky cores are expected to 
be effectively stripped of HSEs after final core segrega- 
tion. However, a surprisingly high abundance of HSEs 
has been inferred on the terrestrial planets, which sug- 
gests continued planetesimal accretio n onto the Earth, 
Moon and Mars after core formati on (| Warren et al.l l 1999: 
iWalker et al.1l200l I Walkerl l2009h . This implies for the 
Earth-Moon system that a significant amount of mass 
was accreted after the Moon- forming impact. The geo- 
chemical evidence for the late veneer and the estimated 
amounts of material added to the Earth, Moon and Mars 
after the giant impact phase are summarized and dis- 
cussed in detail in section 3. 

Independently, there is evidence from planet for- 
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mation models that argues for a population of left- 
over planetesimals in the terrestrial planet region (see 
section [2] for details). This evidence comes from 
planet format ion models that examine the onset of gi- 
ant i mpacts (Ke nvon fc Bromlevl 120061 : iFord fc Chiangl 
l2QQ7h . fro m simulations o f coll isions between pro- 
toplanets (jBenz fc Asphaugl Il999h and from the re- 
quirement to relax the high eccentricities and in- 
clinations of ter r estrial planets after giant impacts 
(iChambersl 120011: lO^Brienet all 120061 : iRavmondetldl 
120061 : iSchlichting fc SarH I2007D . These works argue for 
the existence of a planetesimal population in the terres- 
trial planet region that still contains a few to 10% of the 
mass of the Earth following the epoch of giant impacts. 

In this paper, we investigate the accretion of the left- 
over planetesimals after the end of giant impacts in the 
inner solar system and show that we can account for the 
relative amounts of late veneer accreted by the Earth, 
Moon and Mars provided that most of it was delivered 
by small planetesimals. In section 2 we first estimate 
the mass left in planetesimals after giant impacts from 
planet formation models. We summarize the geochemi- 
cal evidence for a late veneer and estimate the mass of 
late veneer added to Earth, Moon and Mars in section 

3. We calculate the accretion cross section for the Earth, 
Moon and Mars for a range of planetesimals velocities 
and compare our results with the inferred late veneer for 
these bodies from measured HSE abundances in section 

4. We calculate the typical planetesimal sizes that de- 
livered the late veneer in section 5 and estimate their 
accretion timescale in section 6. Discussion and conclu- 
sions follow in section 7. 

2. LEFTOVER PLANETESIMALS AFTER THE END OF 
GIANT-IMPACTS: EVIDENCE FROM PLANET 
FORMATION MODELS 

Orbit crossing and giant impacts begin once mutual 
stirring of the protoplanets can no longer be efficiently 
damped by small planetesimals. Order of magnitude es- 
timates that balance the stirring rates of the protoplanets 
with the damping rates due to dynamical friction, which 
is generated by small planetesimals, find that orbit cross- 
ing sets in when a ~ E, where E and a correspond to the 
mass surface densit y in protoplanets and sm all planetes- 
imals, respectively (jGoldreich et al.ll20 04a). This result 
has been confirmed by numerical simulations studying 
the onset of this instability in the terrestrial planet region 
(jKenvon fc Bromlevl l2006h . Therefore about 50% of the 
total mass still resides in small planetesimals when giant 
impacts begin. During giant impacts planetesimal accre- 
tion continues and additional 'new' planetesimals may 
be produced as byproducts of giant impacts. Therefore, 
due to the fact that giant impacts set in when a ~ E and 
the possible production of 'new' planetesimals in giant 
impacts themselves, a significant population of planetes- 
imals is expected to still be present after the end of giant 
impacts. 

In addition, leftover planetesimals provide a way to re- 
lax the high eccentricities and inclinations of terrestrial 
planets after giant impacts. N-body simulations predict 
eccentricities and inclinations for terrestrial planets af- 
ter giant impacts that are significantly larger than the 
time average d values of the terrestrial p l anets in our 
solar system (Chamb ers fc Wetherilll 119981 : lAgnor et al.l 



1999; Chambers 2001). A population of small planetes- 
imals could provide the dynamical friction that would 
be needed to damp the eccentricities and inclinations 
after giant impacts to the observed values of the ter- 
restrial planets. It has already been shown in direct 
N-body integrations that including a population of less 
massive planetesimals, in addition to the massive plan- 
etary emb ryos, decreases the final eccentricities and in- 
clinations (IChambers fc Wetherilll 1 19981: IChambersl l2001l: 
IQ'Brien et al.ll2006l : iRavmondet al.ll2006D . However, due 
to computational limitations, none of these works were 
able to include planetesimals small enough such that 
their collective interactions could be accurately described 
by dynamical friction. In the limiting case where the ter- 
restrial planets are embedded in a large number of small 
planetesimals, the damping of the velocity dispersion (i.e. 
the damping of the eccentricity and inclination) by dy- 
namical friction only depends on the total mass surface 
density of the planetesimals and is independent of the 
mass of the individual bodies. Since this limit has not 
been reached in direct N-body simulations, such works 
underestimate the strength of dynamical friction for a 
given mass surface density of planetesimals and hence 
overestimate the mass needed in planetesimals to damp 
the eccentricities and inclinations. Further, the veloc- 
ity dispersion of the small planetesimals, u, is likely to 
be damped by mutual planetesimal collisions, which in 
turn may lead to more effective dynamical friction being 
exerted on the terrestrial planets. This is because the 
strength of dynamical friction depends on the relative ve- 
locity between the planetesimals and the terrestrial plan- 
ets, v re i. Collisional damping of the planetesimal velocity 
dispersion has not been modeled by direct N-body simu- 
lations, because including of the order of ~ 10, 000 small 
planetesimals is still not feasible computationally. 

The minimum mass in small leftover planetesimals 
needed to damp the eccentricities and inclinations of the 
terrestrial planets to their current values can be esti- 
mated by comparing the eccentricity damping timescale 
to the accretion timescale of the leftover planetesimals. 
The damping timescale, tdamp, due to dynamical friction 
generated by leftover planetesimals is given by 



dt pR 

tdamp ~ -V— ~ — 



V 
Vesc 



(i) 



where p is the mean density, R the radius, v esc is the es- 
cape velocity of the terrestrial planets and v their velocity 
dispersion (jGoldreich et al.l l2004bh : ft = ^GM Q /a 3 is 
the Keplerian angular frequency around the Sun, where 
Mq is the mass of the Sun and a the semi-major axis. 
The damping timescale needs to be shorter than the time 
required for the remaining planetesimals to be accreted 
by the terrestrial planets. Writing the planetesimal ac- 
cretion timescale as 



dt 
da 



pR 
Eft 



v 

Vesc 



and requiring that t acc > t darnp yields 



cr > E 



(2) 



(3) 



where E corresponds here to the mass surface density of 
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the terrestrial planets. Hence we see from Equation (|3j) 
that we can place a lower limit on the mass required in 
small planetesimals as long as we know the velocity dis- 
persion of the terrestrial planets, v, at the end of giant 
impacts. We assumed in Equations (l)-(3) that v > u, 
where u is the planetesimal velocity dispersion, because 
small planetesimals are likely to damp their velocities by 
mutual collisions ensuring that u <C v. The relative ve- 
locity between the planetesimals and terrestrial planets 
will therefore be dominated by the terrestrial planet's 
velocity dispersion. We can estimate a minimum value 
for v by requiring that the velocity dispersion must have 
been at least large enough for neighboring planetary em- 
bryos that are undergoing giant impacts to cross their 
orbits. N-body simulations find a typical s pacing of plan- 
etary embryos of ~ 8Rh (|Lissaue rT [l993l ). where R H is 
the Hill radius defined as R H = o(M/3M ) 1 / 3 where 
mass of the planetary embryo. This 
suggests that v ~ AVtR H ~ 2.0(R/a) 1 / 2 (M Q /M) 1 / 6 v esc , 
which evaluates to v ~ 0.1v esc for Earth-like terrestrial 
planets at 1 AU. Substituting for v in Equation (J3j) we 
find 



<7>0.01£. 



(4) 



This order of magnitude estimate suggests that at least 
1% of the total mass needs to still reside in small plan- 
etesimals after the giant impact phase in order to damp 
the terrestrial planets' eccentricities and inclinations to 
their observed values. For comparison, N-body simu- 
lations of the giant impact phase, that include smaller 
planetesimals but that are not in the limit in which the 
collective interactions of the small bodies with the ter- 
restrial planets can be accurately described by dynami- 
cal friction, predict eccentr icities of ~ 0-1 f or terrestrial 
planets after giant impacts ([Chamb ers 2001). These sim- 
ulations ther efore suggest v ~ 0.27 v esc and hence that 
a > 0.07S (jSchlichting fc Sarill2QQ7h . Since Equation (gj) 
yields only a lower limit on a and since it is only an order 
of magnitude estimate, we will assume throughout the 
rest of the paper that at least 1% of the total mass needs 
to reside in small planetesimals after the giant impact 
phase in order to damp the terrestrial planets' eccentric- 
ities to their observed values. In addition, we note here 
that the value derived in Equation ((4]) could be reduced 
by a factor of ~ (aQ/Qv) 4 if gaps form in the planetes- 
imal disk (see section 6 for details). This is due to the 
fact that the presence of gaps increases the planetesimal 
accretion timescale, while not significantly altering the 
damping timescale. 

3. GEOCHEMICAL EVIDENCE FOR PLANETESIMAL 
ACCRETION AFTER THE GIANT IMPACT PHASE 

3.1. Earth 

HSEs are comprised of Re, Au and the six platinum- 
group elements Os, Ir, Ru, Pt, Rh and Pd. These ele- 
ments have very high metal-silicate partition coefficients, 
which suggests that the silicate portions of rocky bodies 
with metallic cores should have been stripped of HSEs at 
the end of core formation. Yet, the relative abundances 
of these elements in Earth's mantle are broadly similar to 
chondrites. Absolute concentrations of Ir and Os in the 
Earth's upper mantle are estimated to be ~ 3 — 4ng/g 
(Walker 2009). Although this is more than 100 times 



lower than the concentrations of these elements found in 
chondrites, which typically range from ~ 400 — 800ng/g, 
these concentrations are higher than low pressure part i- 
tion coefficients predict (Horan et al. 2003; Walker 2009). 
Therefore, if the entire mantle harbors HSE abundances 
similar to the estimate for Earth's upper mantle, then 
this suggests, for an Earth mantle mass of ~ 4.0 x 10 27 g, 
that about 1.5 - 4.0 x 10 25 g, or 0.3 - 0.7% M e , of chon- 
dritic material was added to the Earth by late accretion. 

3.2. Mars 

The abundance of HSEs estimated for the Martian 
mantle is ro ughly similar t o that of Earth's prim itive up- 
per mantle (Walker 20091). I Warren et al.l ([1999) estimate 
the primitive mantle abundance of Re, Os, Ir and Au for 
Mars and find that the Martian HSE abundances likely 
range from 0.34-0.66 of the terrestrial values. Assuming 
that the Martian mantle has a mass of 5.1 x 10 26 g, we 
estimate a mass accretion ratio for Earth/Mars of 12-23, 
i.e. the Earth accreted 12 to 23 times more mass as late 
veneer than Mars. In this estimate we assumed that the 
mass accreted by Mars had an HSE composition similar 
to that accreted by the Earth. 

3.3. Moon 

The total amount of late veneer added to the Earth 
can be reasonably well estimated from the HSE abun- 
dances in its mantle but the situation for the Moon may 
be somewhat different, because of the rapid formation of 
a permanent lunar crust by ~ 100 Myr after the Moon- 
forming impact, which contra sts with the delayed devel- 
opment of terrestrial cratons ( Carlso n fc Lugmairl 1988; 
ISnvder et all 120001 : iWalker et al.l 120041 ). Consequently, 
the abundances of HSEs in the lunar mantle probe the 
material accreted before the isolation of the lunar mantle 
by the crust, whereas the HSE abundances in the lunar 
crust provide constraints on the accretion of material af- 
ter crust formation. 

Estimates of the HSE content in the lunar mantle 
have varied consid erably from am ounts similar to that 
of Earth's mantl e (iRingwoocfl 119921) , to amounts about 



20 times lower (IWarren et al 



IDav et all 125571 : [Walker 2(M). However recent works 



1989t IWalker et aLll200i 



generally tend to favor the lower end of thi s range 
(Wal ker et al1l200l IDav et al.ll2QQ7h . Using the |Wker| 
(2009) estimate that the lunar mantle has a factor of 20 
lower HSE concentrations than the terrestrial mantle and 
assuming a lunar mantle mass of 7 x 10 25 g, we find that 
1.3 — 3.5 x 10 22 g of late veneer was added to the lunar 
mantle, yielding an Earth/Moon mass accretion ratio of 
about 1100. However, this estimate ignores an important 
additional HSE repository: the crust and the upper part 
o f the lunar lith o spher e . 

IWalker et al.l (|2004j ) suggested that a significant 
amount of late veneer may have been deposited into the 
lunar crust rather than mantle, implying that the ma- 
jority of the late veneer was accreted by the Moon af- 
ter the formation of a permanent lunar crust. Impact 
melt breccias and bulk regolith sam ples have Os and Ir 
concentrations aver aging 5-15 ng/g ([Morgan et al.lll976t 
iNorman et al.ll2002[ ). Assuming a lunar crustal mass of 
5 x 10 24 g and that the late veneer was delivered by bod- 
ies with Os and Ir concentrations similar to chondrites, 
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which typically range from ~ 400 — 800ng/g, and that 
Os and Ir concentrations are roughly uniform throughout 
the crust, no more than 0.3 — 1.9 x 10 23 g of lat e veneer 
was added to the lunar crust (jWalker et al.ll2004D 0. This 
yields an Earth/lunar crust mass accretion ratio ranging 
from 200 to 700. In this estimate we assumed that the 
Earth and Moon accreted chondritic material with the 
same Os and Ir abundances such that the actual abun- 
dance value cancels in the relative mass accretion ratio. 
The above estimate for the mass added by late accre- 
tion to the Moon assumed, probably somewhat unrealis- 
tically, uniform Os and Ir concentrations throughout the 
lunar crust; below we derive an independent estimate 
of the late veneer added to the Moon by examining the 
eject a layer and the HSE abundance gradient within it. 

The top of the lunar lithosphere is a layer of ejecta 
accumulated from countless impacts. This layer is 
commonly termed a megaregolith, but confusion arises 
because the term megarego lith implies loose debris. 
Pressure-sensitive sintering (jWarrenll2011f ) has probably 
markedly increased cohesion within the deeper portion 
of the ejecta layer. For the loose-debris subvolume of 
the ejecta, a thickness of roughly 2.5 km has been in- 
ferred based on rad ar constraints on blockiness of ejecta 
from large craters (jThompson et al.l I1979L l2QQ9h . But 
the e j ect a- v olume model of W arrenl (j2011f ) combined with 
iFrevI ([2011j Vs inventory of 90 likely lunar basins with di- 
ameters > 300 km suggests a global mean ejecta accu- 
mula t ion of ^5.8 km. Adding the Procellarum basin to 
IFrevI (|2011l ) inventory would increase the accumulation 
by ~ 3.5 km. In Frey's judgment Procellarum shows no 
topographic basin, but if not a single impact, this giant 
region of low elevation and thin crust may reflect a cluster 
of a few ex tremely an cient (degraded), large basins not 
included in IFrevI (|2011[ ). Thus, a reasonable compromise 
basis for estimating the global mean ejecta accumula- 
tion is to add a large fraction, say 1/3, of a Procellarum 
equivalence of thickness, i.e., ~ 1.2 km. This leads to a 
final estimate that the 'known' global ejecta accumula- 
tion thickness is ~ 7 km. However, being based on only 
observable basins > 300 km in diameter, this estimate is 
likely low by a significant factor, of at least order 2. 

The global ejecta layer consists mostly of jumbled tar- 
get (lunar) matter, with a much smaller proportion of im- 
pactor /chondritic matter. For estimating the bulk com- 
position of the ejecta layer, the most useful samples are 
from highly immature regolith. Immature regolith has 
been thoroughly churned and mixed, but not (at least, 
not for long) at the very surface of the Moon, and thus is 
free or nearly free of microme teorite component and as- 
socia ted enrichment in HSE (|McKav et al.lll99it IWarrenl 
2004). For the most commonly measured of the HSE, Ir, 
the average co mposition of all highland regolith samples 
is ~ 12 ng/g (Haskin & Warren 1991). For immature 
highland regolith, with no micrometeorite component, 
the average is probably more like 8 — 10 ng/g. 

The ejecta layer may be only a fraction of the total 
upper-lithosphere component of the Moon's late veneer. 

5 Note, IWalker et all (|200l ) quote a slightly narrower range of 
4 — 8 x 10 22 g for the late veneer, because he assumed a slightly 
narrower range of possible chondritic concentrations of Os and Ir 
than we have used here and assumed a lunar crust al mass of 3.7 x 



Settling of metals in basin-scale impacts has probably 
produced local HSE concentrations deep within and even 
below the crust. An iron-meteoritic or ordinary chon- 
dritic impactor would contain metal as a major mineral. 
In large events, where the central, unejected mass of im- 
pact melt is slow to cool and solidify, the dense metal 
component must tend to settle to the very bottom of 
the impact melt volume. A known example of metal 
that settled out of a lunar i mpact melt is the m ostly 
metallic 4.4 gram rock 14286 (jAlbrecht et al.lll995f ). As- 
suming equilibration occurs, only a tiny proportion of 
metal would suffice to efficiently scavenge the HSE out 
of a silicate melt-metal system. Diffusion within metal 
is very rapid. The limiting factor, for the efficiency of 
HSE scavenging, may be a tendency for the metal com- 
ponents to be so extremely fine grained that they fail to 
settle. However, in the largest events the central mass 
of unejected impact melt is so slow to cool that even its 
silicate crystall ization is believed t o entail gravitationa l 
differentiation (jWarren et all 119961 : llvanov et al.ll201Q[ ). 
The sunken metal probably ended up mostly near the 
bottom of the 'sheet' of central, unejected impact melt, 
at a depth equivalent to roughly 1/10 the diameter of 
the transient crater (jWarren et al.lll99"6h : i.e., in general, 
roughly 1/20 the final basin diameter. In other words, 
the depth at which the metal components predominantly 
settled was probably of order 50 km. It would take a sub- 
sequent basin half as large as the original transient crater, 
if centered at precisely the same point, to begin to exca- 
vate the base of the melt 'sheet'; i.e., the settled metal. 
Thus, HSE concentrations found in the megaregolith at 
the present surface of the Moon may under-represent, 
possibly by a large factor, the total amount of HSE-rich 
matter accreted as late veneer. 

Although large basin-scale impacts likely played a cru- 
cial role in creating and churning the lunar ejecta layer, 
they may not have contributed significantly to the over- 
all lunar HSE budget of the ejecta layer. This is because, 
unlike the small and low velocity impactors that collided 
with the Moon early on and which are the focus of this 
paper, the largest impactors likely collided late in the 
lunar history and had large impact velocities, such that 
only a fraction of the total impact mass was actually 
accreted by the M oon. For example, the Nice model 
(jGomes et al.l l2QQ5h suggests the average Late Heavy 
Bombardment velocity was 21-25 km/s. Under such 
conditions, at the most common impact angles, only a 
fraction of the impactor will actually be accreted by the 
Moon, whereas for the Earth, with its much higher es- 
cape velocity, the fraction of impactor matter that fails to 
accrete is comparatively negligible. Modeling constraints 
suggest that for rocky impactors the lunar accretion ef- 
ficiency, integra ted over all impact angles, i s 0.32-0.16 
for 2 1-25 km/s (jArtemieva fc ShuvalovN2008l : lOng et al.l 

We therefore model the total impact mass that has 
collided with the lunar crust, what may be thought of 
as the equivalent veneer mass, as a combination of two 
main components: the ejecta layer itself, with ~ 10 ng/g 
Ir abundance, and an unobservable, cryptic compo- 
nent. This cryptic component represents the assumed se- 
questered component of metal that settled out as a result 
of larger impacts and, in addition, also accounts for never 
accreted impactor matter that may have been delivered 
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during the Late Heavy Bombardment. The results in 
terms of the Earth to Moon (equivalent) accretion ratio, 
accounting for both the lunar mantle and crust compo- 
nents, are shown in Figure [TJ where the various curves 
correspond to different assumptions for the proportion 
of the cryptic HSE component as a fraction of the total 
lithospheric HSE component. Even if we assume that 
the cryptic component is negligible and that no eject a 
accumulation occurred beyond ~ 7 km from 'known' 
basins with diameters > 300 km, the Earth/Moon ve- 
neer mass ratio is < 700. This result by itself implies that 
the total lithospheric HSE component is important, and 
probably larger than the mantle component, unless the 
sunken-sequestered component is very small. More real- 
istically, the sequestered HSE fraction is probably ~ 0.5, 
but anything between 1/5 and 2/3 seems almost equally 
plausible, and the total eject a accumulation is likely at 
least 2 times the 'known' ejecta volume. The extrap- 
olated Earth/Moon veneer mass ratio is then roughly 
300, assuming 10 ng/g Ir in the ejecta layer. This ra- 
tio would be even lower, if inefficient accretion during 
the Late Heavy Bombardment contributed significantly 
to the cryptic component. The huge uncertainties in the 
sequestered HSE fraction and the contribution of ineffi- 
cient accretion to the cryptic component are permissive 
of the Earth/Moon veneer mass ratio being conceivably 
as low as ~ 150 or as high as ~ 700. 

The implied total masses of chondritic-debris addi- 
tions associated with our various models are perhaps 
most easily comprehended when expressed in terms of 
equivalent proportion of chondritic matter within the lu- 
nar crust (the 'equivalent' qualifier is needed because we 
assume that a major fraction of the HSE actually be- 
came sequestered at the bottoms of the deepest impact 
melt pools near the base of the crust). Assuming the 
mass of the crust is 5 x 10 24 g, models that suggest an 
Earth/Moon accretion ratios of 600, 300 and 150, imply 
equivalent proportions of chondritic matter within the 
crust of about 0.34%, 1.0% and 2.4% by weight, respec- 
tively 0. These proportions represent additions to the 
1.7 x 10 22 g of chondritic matter inferred to be present 
in the mantle source region of the lunar (mare) basalts. 

In summary, we find that the Earth/Moon mass accre- 
tion ratio likely ranges from ~ 150 to ~ 700. In addition, 
we note that comparison of the mass accretion estimates 
for the lunar crust and mantle suggests that most of the 
late veneer was deposited into the lunar crust rather than 
mantle, which indicates that most of the late veneer was 
accreted by the Moon after the formation of a permanent 
lunar crust (jWalker et al.ll2004fl . i.e. about 100 Myrs af- 
ter the Moon- forming impact. 

4. PLANETESIMAL ACCRETION WITH GRAVITATIONAL 
FOCUSING 

4.1. The Accretion of Leftover Planetesimals by the 
Earth and Mars 

If the relative velocity between the planetesimal and 
the terrestrial planets, v re i, is less than the escape veloc- 
ity from the terrestrial planets then the accretion cross 
sections of the terrestrial planets are enhanced above 
their geometric values by gravitational focusing. The 

6 For this illustrative calculation, we assume that all the material 
delivered to the Earth and Moon was accreted. 
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Figure 1. Earth/Moon (equivalent) accretion ratio, shown as a 
function of assumed mean global thickness of the ejecta layer. The 
various curves are labeled to indicate different proportion of the 
cryptic HSE component as a fraction of the total lithospheric HSE 
component. The cryptic component represents the assumed se- 
questered HSE component and, in addition, accounts for never- 
accreted impactor matter that may have been delivered during the 
Late Heavy Bombardment. The ejecta layer is assumed to contain 
10 ng/g Ir. 

gravitationally enhanced cross section is given by 

2" 



A = ttR 



1 



/ Vesc 
\Vrel 



(5) 



Using Equation we can write the accretion cross sec- 
tion ratio of Earth/Mars as 



(1 + (v esc (®)/v re i) 2 )Rl 
(1 + (v esc (Mars)/v rel ) 2 )R 2 Ma 



(6) 



where the subscripts and 'Mars' label the quantities 
corresponding to Earth and Mars, respectively. If v re \ <C 
v esci Equation [6] can be simplified to 



pMar 



Rn 



Rm a 



17, 



(7) 



where p e and pMars correspond to the mean density of 
the Earth and Mars, respectively. This implies that the 
accretion ratio between the Earth and Mars that has 
been estimated from the HSE abundances (A^/AMars ~ 
12 — 23, see section 3) agrees very well with the expected 
mass accretion ratio between the Earth and Mars, if the 
accretion cross sections of the Earth and Mars were sig- 
nificantly enhanced by gravitational focusing. For com- 
parison, if v re i had been larger than the escape velocities 
of Earth and Mars such that gravitational focusing be- 
comes irrelevant, then the ratio of the accretion cross sec- 
tion Earth/Mars is simply given by (R^/RMars) 2 ~ 4. 
This value is significantly lower than the mass ratio of 
the late veneer that has been estimated from the terres- 
trial and Martian abundances of HSEs (see Figure [2]). 
We therefore conclude that we can account for the rel- 
ative amounts of late veneer accreted by the Earth and 
Mars, if it was delivered concurrently by a population 
of planetesimals with a velocity dispersion small enough 
such that v re i <C v esc . 
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Figure 2. Ratio of the mass accretion rates for the Earth and 
Mars, /AMars-, as a function of v re \ in units of Earth's es- 
cape velocity, v eS c((B)- The dotted line represents the limit with- 
out any gravitational focusing (v re i ^> u e sc(©)), in which case 
A&/AMars — RMars) 2 ~ 4. The dashed line corresponds 

to the limit with strong gravitational focusing (v re i <C u eS c(©))) 
in which case A e /A M ar S = 0© / 'PMars)(R® / '^Mars) 4 ~ 17. The 
yellow shaded region shows the range of Earth/Mars mass accre- 
tion ratios that are consistent with the inferred late veneer for the 
Earth and Mars. 

4.2. The Accretion of Leftover Planetesimals by the 
Earth and the Moon 

We can extend the above argument to the Earth and 
the Moon to get a rough estimate for the Earth/Moon ac- 
cretion ratio. However, using the expression for gravita- 
tional focusing from Equation [6] is strictly speaking only 
valid for isolated bodies and is therefore only a rough 
approximation for the Earth/Moon accretion ratio. Es- 
timating the Earth/Moon accretion ratio from Equation 
[6] we have 



^Moon 



pMo 



Ra 



Rmo 



300. 



(8) 



This implies that in the limit in which we can treat the 
Earth and Moon as isolated bodies, i.e., for very large 
Earth-Moon separations, the Earth should have accreted 
about 300 times more mass as late veneer compared to 
t he Moon (see Figure 03). 

iBandermann fc Singer! (|1973l ) derived analytically the 
ratio of the Earth/Moon accretion cross section. As- 
suming an isotropic planetesimal velocity distribution far 
from the Earth and neglecting Earth's shadow they find 



An 



A 



> 



\ + {v rel /v esc (®)f 



Moon 



l(R®/a E -M) + 0.045 + (Vrel/Vesc(®)) 2 \Rmoot, 



Ra 



(9) 

where ole-m is the Earth-Moon separation. Equation (j9j) 
is a lower limit to A^/Amoou because the effect of Earth's 
shadow was neglected in deriving the accretion cross sec- 
tion for the Moon, i.e., it neglects lunar impactors that 
would have collided with Earth first. We confirmed this 
analytic result for the Earth/Moon accretion cross sec- 
tion by direct numerical integrations of planetesimal tra- 
jectories in the Earth-Moon system (see Figured]). 

Figure [3] shows the Earth/Moon accretion ratio as a 
function of the Earth-Moon separation for v re \ <C v esc . 
As expected the Earth-Moon accretion cross section ap- 




500 



Figure 3. Ratio of the mass accretion rates for the Earth 
and Moon, A^/Amoou, as a function of the Earth-Moon dis- 
tance in Earth radii, cle-m / R@ with strong gravitational focus- 
ing (i.e., v re \ <^ v e sc) : The solid line is the analytic result from 
Bandermann & Singer (1973) given in Equation (0. The dashed 
line represents our estimate for the ratio of the Earth/Moon accre- 
tion cross section from Equation [8] which is valid for large Earth- 
Moon separations when the Earth and the Moon can be well ap- 
proximated as isolated bodies. The yellow shaded region shows 
the range of Earth/Moon mass accretion ratios that are consistent 
with the inferred late veneer added to the Earth and Moon based 
on HSE observations. 



proaches the ratio for isolated bodies at large Earth- 
Moon separations. At the current Earth-Moon separa- 
tion, which corresponds to about 60 Re, the accretion 
ratio is ~ 200. The Earth/Moon accretion ratio increases 
with increasing Earth-Moon separation approaching the 
limit derived in Equation [8] for isolated bodies, this be- 
havior remains unchanged as long as gravitational focus- 
ing plays a significant role in enhancing the accretion 
cross section above the geometric value (see Figure H]). 
This result may seem surprising at first, because Earth's 
gravitational field accelerates the incoming planetesimals 
such that the Moon's gravitational focusing is reduced. 
However, Earth's gravitational field also focuses the in- 
coming planetesimals such that the Moon intercepts a 
larger number of planetesimals than it would have other- 
wise. As a result, the Earth/Moon accretion cross section 
decreases with decreasing Earth-Moon separation, which 
implies that it was lower in the past before the Moon 
evolved tidally outward to its current location. The tidal 
evolution timescale for the Moon to evolve from an initial 
separation of a few times Earth's radius to a separation 
%e-m > R® is 

a E _ M Y(GMv_yV\ L3xlQ4 



tidal 



2 Q Mg 
39 k M Mo 



Ra 



\a 3 



E-M 



(10) 

where k and Q are the tidal dissipation function and 
the tidal Love number of the Earth, respectively. From 
Equation [10] we see that the initial tidal evolution was 
very fast such that the Earth-Moon system did not spend 
a significant amount of time at small Earth-Moon sep- 
arations and evolved to cle—m ^ 40i?e within about 
110 Myrs of the Moon forming impact. This estimates 
assumes that Q ~ 12 and that its value did not signif- 
icantly change throughout the tidal evolution. In sec- 
tion 6 we calculate a planetesimal accretion timescale 



Q/12 
/c/0.299 
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of ~ 170 Myr, which suggests that most of the plan- 
etesimals were accreted after the Earth-Moon system 
evolved to separations > 40i? e . We therefore conclude 
that the relative amounts of late veneer added to the 
Moon and Earth, as inferred from their HSE abundances, 
are consistent with the accretion of small planetesimals 
with a velocity dispersion u < v re i ~ O.li> esc (0). In 
contrast, if v re i ^> v esci i.e. in the limit without any 
gravitational focusing, the Earth/Moon accretion ratio 
is independent of the Earth-Moon separation and is sim- 
ply given by the ratio of the geometric cross sections 
(R^/Rmoou) 2 ~ 14 (see Figure Hj), which is inconsistent 
with the Earth/Moon mass accretion ratio inferred from 
HSEs of 150-700. 




WVesc(©)=2 



10 20 30 40 50 60 70 

&e-mI ft© 

Figure 4. Ratio of the mass accretion rates for the Earth and 
Moon, A© /A M com as a function of the Earth-Moon distance in 
Earth radii, cle —m/R^ for various relativ e velocities, v re i, from 
Equation J5J (Bandermann & Singer 1973). For comparison, we 
also show the results from our numerical simulations (shown as 
points) in which we integrate the trajectories of the planetesimals in 
the Earth-Moon system directly. The yellow shaded region shows 
the range of Earth/Moon mass accretion ratios that are consistent 
with the inferred terrestrial and lunar late veneer based on HSE 
observations. 



damping rates yields 

^Ps \VescJ \VescJ 

Evaluating Equation ([T2]) for cr/E ~ 1% (see discussion 
in section 2), substituting for v re \ ~ v ~ 0.1v esc and 
using u ~ v, we have 

s ~ 10 m (13) 

where we assumed p s ~ 3 g/cm 3 and R ~ R® and 
p ~ p0. This implies that u < v ~ 0.1v esc (©) as long as 
the typical planetesimals, that damped the eccentricities 
and inclinations of the terrestrial planets and that deliv- 
ered the late veneer, were smaller than about 10 meters 
in size. We assumed when evaluating Equation ([T3]) that 
u ~ v. The actual planetesimal sizes could therefore be 
smaller than estimated in Equation ([13]) . which would 
imply that u < v. The corresponding optical depth of 
such a planetesimal population is r ~ cr/p s s, which im- 
pli es r > 2.5 x 10~ 5 : We also note here that recent work 
by IWeidenschillingl (|2Qllh showed that the size distribu- 
tion of the asteroid belt can be reproduced by coagula- 
tion from an initial population of planetesimals as long as 
they have sizes < 100 m, providing independent support 
for a planetesimal population in the inner solar system 
that was smaller than about 100 meters in size. 

If typical planetesimal sizes would have exceeded about 
10 meters, then u > v such that v re \ ~ u > 0.1v esc , which 
would imply weaker or no gravitational focusing, mak- 
ing it hard to reconcile the resulting Earth/Moon and 
Earth/Mars accretion ratios with the relative quantities 
of late veneer delivered to these bodies. Furthermore, if 
u > v ~ 0.1v esc than this would imply that more mass, 
than we estimated in section 2, must have been resid- 
ing in small planetesimals at the end of giant impacts in 
order to damp the eccentricities and inclinations of the 
terrestrial planets. 



5. PLANETESIMAL SIZES 

We have shown in section 4 that the relative amounts 
of late veneer delivered to Earth, Moon and Mars can be 
explained by the accretion of small planetesimals after 
giant impacts provided that gravitational focusing played 
a significant role in increasing the accretion cross section 
above the geometric value. Since gravitational focusing 
only acts when v re \ <C v esci we conclude that the velocity 
dispersion of the planetesimals was comparable to or less 
than the velocity dispersion of the terrestrial planets, i.e. 
u < such that v re i ~ v. Further, we can use our 
estimate from section 2 that v ~ 0.1v esc to place an 
upper limit on the typical planetesimal size as follows. 

The velocity dispersion of the planetesimals is stirred 
gravitationally by the terrestrial planets and damped by 
mutual planetesimal collisions such that 

1 du YiVt ( v esc \ fv esc 

U dt pR \V re l J V U J p s S 

where s and p s are the radius and density of the plan- 
etesimals and u < v re i. Balancing the stirring and the 



6. GAP FORMATION AND ACCRETION TIMESCALES 

Gaps were not important before the end of the gi- 
ant impact phase, because the radial separation of pro- 
toplanets was only a few times larger than the widths 
of their Hill radii. But at the end of giant impacts, 
when the terrestrial planets achieved large-scale orbital 
stability, their orbital separation was much larger than 
their Hill radii and gaps lik ely formed around their or- 
bits (|Goldreich et al ll2004af ). Gaps increase the clean up 
timescale of the leftover planetesimals, because accretion 
onto the protoplanet s can now only proceed from the 
gap edges. Following llGoldreich et all (|2004a| ). the rate 
at which a terrestrial planet accretes small planetesimals 
from gap edges located a distance x from the planets' 
semi-major axis is 

M dt pR \Aa J \vrei J 

where Aa is the distance between neighboring plan- 
ets. The gap surface density of the small bodies fol- 
lows a power law such that the mass surface density at 
the gap edges is given by a = <Jo(2x/Aa) 4 . Writing 
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tacc = -^o^ ~ 2naoaAa(dM/dt)- 1 yields: 




where E = M/(27raAa). The spacing between the ter- 
restrial planets is Aa ~ a/3 and, since v > u, x will 
be roughly given by the radial excursions of the planet 
which is ~ ae. This yields a clean up timescale of 

t acc = ~ ^ f^V^Y - 170 Myrs 

Ct<7 SO VVesc/ \6e/ 

(16) 

when evaluated at 1 AU for R = R®, M = M e and 
v re i ~ eaQ ~ 0.1v esc . This implies that most of the plan- 
etesimals were accreted after the formation of a perma- 
nent lunar crust. This result is consistent with the fact 
that the majority of the late veneer seems to reside in the 
lunar crust rather than mantle (see section 3 for details). 
Because the planetesimal accretion timescale is long com- 
pared to the tidal evolution timescale of the Earth-Moon 
system (see Equation (10) in section 4), this implies that 
the Earth-Moon separation was already > 40 R® when 
most of the late veneer was delivered to the Earth-Moon 
system. 

7. DISCUSSION AND CONCLUSIONS 

The abundances of HSEs suggest that a total of about 
0.01 M0 was delivered as 'late veneer' by planetesimals 
to the terrestrial planets after the giant impact phase. 
We showed here that small residual planetesimals, with 
radii < 10 m, containing about ~ 1% of the mass of the 
Earth can provide the dynamical friction needed to re- 
lax the terrestrial planets' eccentricities and inclinations 
after giant impacts and simultaneously account for the 
relative and absolute amounts of late veneer added to 
Earth, Moon and Mars. Small planetesimal sizes are 
required to ensure efficient damping of the planetesi- 
mal's velocity dispersion by mutual collisions, which in 
turn ensures sufficiently low relative velocities between 
the terrestrial planets and the planetesimals such that 
the planets' accretion cross sections are significantly en- 
hanced by gravitational focusing above their geometric 
values. Specifically we find that, if v re \ <C v esci gravita- 
tional focusing yields a mass accretion ratio Earth/Mars 
mars )(R@/Rmars) ~ 17, which agrees well with 
the mass accretion ratio inferred from HSEs of 12-23. For 
the Earth-Moon system, we find a mass accretion ratio of 
~ 200, which, as we show in section 3, is consistent with 
estimates of 150-700 derived from HSE abundances that 
include the lunar crust as well as mantle component. Fur- 
ther, we find that the higher abundance of siderophilic 
elements in the lunar crust compared to the lunar mantle 
is consistent with the idea that most of the late veneer 
was delivered by small planetesimals. This is because, 
if the residual planetesimals were indeed small, gaps will 
likely form around the terrestrial planets, which will pro- 
long the planetesimal accretion timescale such that most 
of the late veneer is added to the lunar crust after the 
lunar mantle was isolated by the formation of a perma- 
nent crust. We note here that, although we suggest that 
the majority of the accreted late veneer was delivered by 
small planetesimals, larger planetesimals were certainly 
residing among the small planetesimal population and 



must have played a crucial role the mixing and settling 
o f HSEs in the l unar crust. 

iBottke et al.l ([2010) suggested recently that most of 
the late veneer may have been delivered by a few very 
large planetesimals with the largest terrestrial impactor 
exceeding more than 1000 km in radius. Delivering the 
majority of the late veneer by one or two very large 
bodies may explain stochastically a large mass accretion 
ratio between the Earth and Moon. However, whereas 
small body accretion can account for the relative quanti- 
ties of late veneer inferred from HSE abundances for the 
Earth/Moon and Earth/Mars simultaneously, it would 
remain a coincidence in a stochastic accretion scenario. 
Further, if a small number of 1000 km sized planetes- 
imals were indeed responsible for the late veneer, then 
these planetesimals cannot have damped the eccentrici- 
ties and inclinations of the terrestrial planets after giant 
impacts. This is because, if most of the planetesimal 
mass resided in such large bodies, they would have to 
have a significantly higher velocity dispersion, because 
mutual planetesimal collisions, that damp their veloci- 
ties are significantly less frequent for larger planetesimals 
compared to small ones (see section 5). In this case v re \ 
would be determined by the velocity dispersion of these 
large planetesimals, which in turn implies that signifi- 
cantly more than 1% of the total mass would be required 
in large planetesimals to damp the eccentricities and in- 
clinations of the terrestrial planets (see section 2). This, 
however, would be inconsistent with the < 0.01 of 
chondritic material delivered as late veneer to the Earth, 
Moon and Mars. In principle a population of planetesi- 
mals made of primarily silicates with extremely low HSE 
abundances could have provided the required dynamical 
friction. However, the relative abundances of the dif- 
ferent HSEs in the terrestrial planets and the Moon are 
consistent with chondritic material and hence favor the 
idea that they were delivered by small, undifferentiated 
planetesimals with chondritic composition. 

Finally, as we have shown in section 3, the Earth/Moon 
impact ratio likely falls in the range 150-700, once the 
HSE deposited into the lunar crust and the upper part 
of the lunar lithosphere are accounted for (see also 
Wal ker et al.l (|2004l )). If the ratio of the late veneer ac- 
creted by the Earth and Moon falls at the lower end of 
this range then it is consistent with small body accre- 
tion. However, if it can be conclusively shown that the 
Earth/Moon mass accretion ratio lies at the upper end of 
this range then it cannot be explained by the small body 
accretion discussed here. In this case it could instead be 
either due to a small number o f stochastic events t hat de- 
livered most of the late veneer ( Bo ttke et al.ll2010[ ) or due 
to a smaller retention fraction of the material delivered 
to the Moon compared to Earth as might be expected for 
impact velocities significantly exceeding the lunar escape 
velocity. 



We thank David Jewitt and our two referees, John 
Chambers and Richard Walker, for valuable comments 
that helped to improve this manuscript. For HS sup- 
port for this work was provided by NASA through Hub- 
ble Fellowship Grant # HST-HF-51281.01-A awarded 
by the Space Telescope Science Institute, which is op- 
erated by the Association of Universities for Research 



The Last Stages of Terrestrial Planet Formation 



9 



in Astronomy, Inc., for NASA, under contact NAS 5- 
26555. PW acknowledges support from NASA grant 
NNX09 AE31 G. QZY acknowledges NASA Cosmo- 
chemistry (NNX11AJ51G) and Origins of Solar Systems 
(NNX09AC93G) grants for support. 



REFERENCES 



Agnor, C. B., Canup, R. M., & Levison, H. F. 1999, Icarus, 142, 
219 

Albrecht, A., Herzog, G. F., Klein, J., Middleton, R., Schultz, L., 
Weber, H. W., Kallemeyn, G. W., & Warren, P. H. 1995, in 
Lunar and Planetary Institute Science Conference Abstracts, 13 

Artemieva, N. A. & Shuvalov, V. V. 2008, Solar System Research, 
42, 329 

Bandermann, L. W. & Singer, S. F. 1973, Icarus, 19, 108 

Benz, W. & Asphaug, E. 1999, Icarus, 142, 5 

Bottke, W. F., Walker, R. J., Day, J. M. D., Nesvorny, D., & 

Elkins-Tanton, L. 2010, Science, 330, 1527 
Carlson, R. W. & Lugmair, G. W. 1988, Earth and Planetary 

Science Letters, 90, 119 
Chambers, J. E. 2001, Icarus, 152, 205 
Chambers, J. E. & Wetherill, G. W. 1998, Icarus, 136, 304 
Chiang, E. & Youdin, A. N. 2010, Annual Review of Earth and 

Planetary Sciences, 38, 493 
Day, J. M. D., Pearson, D. G., & Taylor, L. A. 2007, Science, 315, 

217 

Ford, E. B. & Chiang, E. I. 2007, ApJ, 661, 602 

Frey, H. 2011, in Recent Advances and Current Research Issues 

in Lunar Stratigraphy, ed. Ambrose, W. A. and Williams, D. 

A., 53-75 

Goldreich, P., Lithwick, Y., & Sari, R. 2004a, ApJ, 614, 497 
— . 2004b, ARA&A, 42, 549 

Gomes, R., Levison, H. F., Tsiganis, K., & Morbidelli, A. 2005, 

Nature, 435, 466 
Haskin, L. A. & Warren, P. H. 1991, in Lunar Sourcebook, A 

User's Guide to the Moon, ed. Heiken, G., Vaniman, D., and 

French, B. M., 357-474 
Horan, M., Walker, R., Morgan, J. W., & Grossman, J. N. 2003, 

Chemical Geology, 196, 5 
Ida, S. & Makino, J. 1993, Icarus, 106, 210 



Ivanov, B. A., Melosh, H. J., &; Pierazzo, E. 2010, in Large 
Meteorite Impacts and Planetary Evolution IV, ed. Gibson R. 
L. and Reimold, W. U., 29-49 

Lissauer, J. J. 1993, ARA&A, 31, 129 

McKay, D. S., Heiken, G., Basu, A., Blanford, G., Simon, S., 

Reedy, R., French, B. M., & Papike, J. 1991, in Lunar 

Sourcebook, A User's Guide to the Moon, ed. Heiken, G., 

Vaniman, D. and French, B. M., 285-356 
Kenyon, S. J. &; Bromley, B. C. 2006, AJ, 131, 1837 
Morgan, J. W., Gros, J., Takahashi, H., & Hertogen, J. 1976, in 

Lunar and Planetary Science Conference Proceedings, ed. 

R. B. Merrill, 2189-2199 
Norman, M. D., Bennett, V. C, & Ryder, G. 2002, Earth and 

Planetary Science Letters, 202, 217 
O'Brien, D. P., Morbidelli, A., & Levison, H. F. 2006, Icarus, 184, 

39 

Ong, L., Asphaug, E. I., Korycansky, D., & Coker, R. F. 2010, 
Icarus, 207, 278 

Raymond, S. N., Quinn, T., & Lunine, J. I. 2006, Icarus, 183, 265 
Ringwood, A. E. 1992, Earth and Planetary Science Letters, 111, 
537 

Schlichting, H. E. & Sari, R. 2007, ApJ, 658, 593 

Snyder, G. A., Borg, L. E., Nyquist, L. E., & Taylor, L. A. 2000, 

in Origin of the Earth and Moon, ed. Canup, R. M., Righter, 

K., & et al., 361-395 
Thompson, T. W., Campbell, B. A., Ghent, R. R., & Hawke, 

B. R. 2009, Geology, 37, 655 
Thompson, T. W., Roberts, W. J., Hartmann, W. K., Shorthill, 

R. W., &; Zisk, S. H. 1979, Moon and Planets, 21, 319 
Walker, R. J. 2009, Chemie der Erde / Geochemistry, 69, 101 
Walker, R. J., Horan, M. F., Shearer, C. K., & Papike, J. J. 2004, 

Earth and Planetary Science Letters, 224, 399 
Warren, P. H. 2004, in Treatise on Geochemistry, Volume 1, 

Meteorites, Comets, and Planets, ed. Davis, A. M., 559-599 
Warren, P. H. 2011, Meteoritics and Planetary Science, 46, 53 
Warren, P. H., Claeys, P., & Cedillo-Pardo, E. 1996, in The 

Cretaceous-Tertiary Event and Other catastrophes in Earth 

History, ed. Ryder, G., Fastovsky, D. and Gartner, S., 105-124 
Warren, P. H., Jerde, E. A., & Kallemeyn, G. W. 1989, Earth and 

Planetary Science Letters, 91, 245 
Warren, P. H., Kallemeyn, G. W., & Kyte, F. T. 1999, 

Geochim. Cosmochim. Acta, 63, 2105 
Weidenschilling, S. J. 2011, Icarus, 214, 671 
Weidenschilling, S. J., Spaute, D., Davis, D. R., Marzari, F., & 

Ohtsuki, K. 1997, Icarus, 128, 429 



